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Abstract. It is well known that some quantum and statistical fluctuations of a 
quantum field may be recovered by adding suitable stochastic sources to the mean 
field equations derived from the Schwinger-Keldysh (Closed-time-path) effective action. 
In this note we show that this method can be extended to higher correlations and 
higher (n-particle irreducible) effective actions. As an example, we investigate three 
and fourth order correlations by adding stochastic sources to the Schwinger - Dyson 
equations derived from the 2-particle irreducible effective action. This method is a 
simple way to investigate the nonlinear dynamics of quantum fluctuations. 

PACS numbers: 03.70.+k, ll.10.Wx 
1. Introduction 

Quantum fields fluctuate, and quantum fields out of equilibrium show both quantum 
and statistical fluctuations In many problems of interest, the fluctuations are more 
relevant than the mean fields themselves. Problems that come to mind are the generation 
of primordial fluctuations during inflation [2j [3J, SI El [6] , the fluctuations of soft fields 
induced by the interaction with hard quanta EJ 0, [101 EJ EE], and the fluctuations 
of a Bose-Einstein condensate as described by the stochastic Gross-Pitaievskii equation 
[T3| [HI [151 EE EH EH [E] ■ hi suc h a case, one can try to obtain information about the 
average behavior of the fluctuations by deriving equations of motion for the fluctuation- 
fluctuation correlations, or else one may attempt to investigate the space-time unfolding 
of the fluctuations by deriving suitable Langevin-like equations for them. In certain 
cases, the Langevin approach is also an efficient way to derive the required self- 
correlations [20J. 

In the simplest set up, one is dealing with a bosonic field theory. The Heisenberg 
field operators are designed as We use a De Witt notation where a accounts for both 
discrete and space-time indexes; repeated indexes are summed over the discrete ones 
and integrated over the continuous ones. The $^ have expectation values = <f) a - If 
only the mean fields are relevant, we may obtain causal equations of motion from them 
from the Schwinger-Keldysh (or closed-time-path (CTP)) 1-particle irreducible (1PI) 
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effective action (EA). Because the Schwinger-Keldysh approach involves doubling the 
degrees of freedom, an extra discrete index appears, and the fields become $^ and <p A . 
In the simplest representation, A = (i,a), where i = 1,2 shows in which branch of the 
CTP we are. Other representations are also possible. The physical mean field equations, 
notwithstanding, are obtained from the CTP equations by adding the constraints 

<f) la = <J) 2a = <J) a (1) 

If fluctuations are important, one may add noise terms to the physical mean field 
equations of motion. The noise self correlation is also derived from the 1PI CTP EA, 
more precisely from terms that do not contribute to the mean field equations when the 
constraints (pQ) are enforced (see next Section). The resulting theory is still good enough 
to derive exact symmetric expectation values for the product of two quantum fields, that 
is, the Hadamard propagator 

Gf =({$ a H , ® h H })-W<t> h (2) 

and in this sense it is a nontrivial extension of the mean field theory. This basic 
framework has been extensively used in cosmology (see pQ; some influential papers are 
[2H E21 [23]; see also the reviews [21 [25]; for more recent work see [26l M, M, EHJ [30] ) 
and in the theory of Bose-Einstein condensates. The important point of to which extent 
these fluctuations may be considered real is discussed in pQ; for present purposes, it is 
enough to consider the stochastic approach as a shortcut to the actual propagators. 

Sometimes higher correlations are also important. For example, one may want 
to compute the expectation value of the energy-momentum tensor (or the fluctuations 
thereof) in an interacting bosonic field theory -that usually involves three and four 
point functions [3H [32J [331 EH E3, [363 • Density - density correlations in a Bose-Einstein 
condensate are a four point correlation of the fundamental Heisenberg field; these are 
relevant, for example, when the condensate is investigated through Bragg scattering 
|37i [38j [39]. One may need accurate higher order order correlations to enforce important 
Ward-Takahashi or Slavnov- Taylor identities [101 [H]. Or simply one may want to 
compute higher cumulants as a way of accounting for nonlinear effects |42j . In this case, 
one of the most powerful computational tools is to obtain self-consistent Schwinger- 
Dyson equations from variations of higher n-particle irreducible (nPI) EAs where all 
required correlations appear on equal footing as independent variables [HU HU H3, HE] . 

There are both practical and fundamental reasons to wish to extend the stochastic 
approach to field fluctuations to higher correlations as well [HJ EE] . On the practical 
side, getting the full fourth order correlations from a stochastic approach to the 2PI EA 
may be more efficient (or at least more heuristic) than computing the whole 4PI EA. 

On the fundamental side, let us mention the following issue. It is well known that 
the 2PI equations of motion for the propagators lead to the Kadanoff-Baym equations 
for the density of states and one-particle distribution function, and eventually to the 
Boltzmann equation (or similar) in the appropriate limit. However, it is also known that 
the Boltzmann equation is only a mean field approximation to a stochastic equation, 
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where the noise terms, in the near equilibrium case, may be derived from the fluctuation- 
dissipation theorem. Of course, field theory complies with the Kubo-Martin-Schwinger 
theorem, and therefore has the fluctuation-dissipation theorem built in. So the noise 
terms in the stochastic Boltzmann equation must correspond to some elements already 
present in the 2PI EA. The fundamental question is to make those elements explicit. 
This issue was solved by Hu and one of us in ref. [19]. Similar issues appear at every 
order in the Schwinger- Dyson hierarchy. 

The stochastic approach to the 1PI CTP EA exploits special features of this EA 
(see below) and is not readily generalizable to higher EAs. Similarly, the approach of 
Calzetta and Hu in ref. [19] is also ad hoc, in this case for the 2PI EA. A systematic 
framework, which could be applied to any EA and to the symmetry broken or unbroken 
cases alike, would be highly desirable, not least because of the light it sheds on the 
particular approaches devised for the 1PI and 2PI cases. Our aim in this note is to 
develop such an uniform formalism. 

The rest of the paper is organized as follows. 

In the next section we study noise in the 1PI theory. We first show how one can 
associate a CTP 1PI EA efective action to a problem defined in terms of a Langevin 
equation. We then show that the 1PI EA for a quantum field theory problem has, under 
certain approximations, the same structure as the EA arising from a stochastic problem. 
This allows the direct identification of the equivalent stochastic problem to a given field 
theory. As an application, we review the derivation of the Hadamard propagator of the 
full theory from the equivalent stochastic equation. 

In the following Section we review the 2PI EA and two early attempts of a stochastic 
formulation of the propagator dynamics [48, 49J. We show the shortcomings of these 
attempts and how they differ from the proposal in this note. 

Finally we present a systematic approach to building stochastic equivalents for a 
given field theory and apply it to the 1PI and 2PI cases. Although we shall not discuss 
it explicitly, generalization to higher effective actions is straightforward. In the 2PI 
case, we finally obtain the same result as in [19], but without the contrived arguments 
contained in that paper. We show the basic oversight contained in jl9], which obscured 
the simple derivation of the 2PI noise presented here. 

The paper ends with some brief final remarks. 

2. Stochastic approach to the 1PI EA 

The goal of this Section is to provide an heuristic introduction to stochastic equations 
derived from the 1PI EA. For a deeper discussion see [Tl 150]. 

2.1. From Langevin equations to effective actions 

To see why it is natural to translate a problem described in terms of an effective action 
into an equivalent Langevin equation framework, let us first traverse the opposite road, 
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that is, how to associate a generating functional to a problem whose primary description 
is in terms of a stochastic equation of motion. 

We therefore assume we have a string of c-number fields $" obeying a system of 
equations of the form [5T| 1521 153} IMII55] 

D a [$J = -g b F a 6 [$J (3) 

where D a and i 7 ^ are possibly nonlinear functionals and the g& are stochastic Gaussian 
variables with zero mean and self-correlation (qbq c ) = Qbc- To avoid the complexities 
associated to nonlinear Langevin equations, we take this equation to mean that it is 
possible to write = cf) a + (p a , where <p a is a solution of the homogeneous equation 

D a [0] = (4) 

and 

D a , c [0]^ = -g fe F a b [0] (5) 
where here and henceforth a comma denotes a functional derivative 

D ac = ^ (6) 

a ' c 8(p c 1 ; 

If there are no zero modes, this implies that the expectation value {(p a ) = 0, so within 
this approximation we may say that a is the expectation value of 

The generating functional for expectation values of the fields is 

Z [J] = e lW[J] = J D& s DqQ5 ($ a s - a - G* b et q c F£ [0]) e iJ ^° (7) 

where G^ b et is the causal Green function for the operator D a c , namely 

D a , c Gf et = -8 b a (8) 

and Q is the Gaussian probability density functional of the q sources. The generating 
functional obeys 

W [0] = (9) 
SW 



W' a [0] 



(10) 



J=0 



Shifting the <3>" fields by an amount a , namely = a + , we get 

e iW[J] = e iJA d J Dip+ DqQ6 ^* _ Gf et q c F b c [0]) e* Jd < (11) 



This may be rewritten as 

JJ d <f> d 



e lW[J] = J D V+ Dq Q 5 (d«^ + q c F c a [0]) e iJ *< (12) 
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We now introduce an auxiliary field ip_ to exponentiate the delta function 



,iW[J] 



jj d <j> d 



DetG re t 

and perform the functional integral over the sources 

P iJd<t> d 



j Dip + Dip_ Dq Q e ^-( D ^++^ F »])e iJ ^+ 



JW[J] 



DetG ret 



where 



N ab = Q cd F c a [0] F b d [0] 



(13) 



(14) 



(15) 



This generating functional is the particular case, when J + = J and J = 0, of the 
functional 



,iw[j+,j~] = e d P 



DetG 



ret 



j Dip + Dip.e^'^+-^- N ^-A^< +J ^-\ (16) 



Let us write 



W 



J + ,J~ 



JU d + w 



J + ,J- 



(17) 



w [J + , J ] is obtained from a functional Fourier transform of a Gaussian functional, and 
as such it is a quadratic function. It therefore obeys the Euler theorem 



w 



where we have introduced two new background fields 



d 5w 
Sw 



sj: 



(18) 

(19) 
(20) 



The expectation value of the + field in the presence of non zero external sources is 

<i> d + = jj+ = <!> d + <p d + (21) 

We can now compute the effective action as the Legendre transform of the 
generating functional 



r [0 + , 0_] = w r] - { j+< + (/>*_} 

We have the identities 

D„ r l = —jf 



-1 



{J+Vi + J^t} (22) 



(23) 
(24) 



and therefore 



+ ^<l> d -N db( t> b _ 



(25) 
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or else, adding D a [</>] =0 
r \6>, 6 J = 



> + } + -<f> d LN db [0 + ] ^ 



(26) 



This shows how we can associate an effective action to a Langevin type equation. 

We finally mention the self correlation for the (p a fields. Since we have the explicit 
representation ip a = G°^ et q c F^ we get 



(27) 



2.2. From effective actions to Langevin equations 



Let us now investigate a scalar field theory described by the Heisenberg operators <& a H . 
The fields have one and two-particle expectation values 



G 



ah 



(28) 
(29) 



These expectation values cannot be derived either from the euclidean generating 
functional or its analytic continuation to Minkowski space, which generate IN-OUT 
matrix elements instead p] . To find a suitable generating functional, we must consider 
two external sources J A = ( J la , J 2a ) and introduce the Schwinger - Keldysh or closed 
time-path (CTP) 1 particle-irreducible (1PI) generating functional W\pi [Ja] as 



npi 



[J] 



D iW 1PI [J] 



(f 



g ' J a. ^ H 



(T 



(30) 



Where T (T) means (anti) temporal ordering. We may now derive the expectation value 

0J a J1=J2 =0 

More generally, we may consider the expectation value of the Heisenberg operator driven 
by an external source J a . This is 



5W 



ipi 



(32) 



JW 2 =J 

Even more generally, we may consider this as the theory of a field doublet $ A = 
($ la , $ 2a ) coupled to sources Ja = (Ji a = J la , -ha = —J 2a ) (observe that lowering or 
raising a 2 index involves a sign change). In this theory we have two backgrounds fields 

lA 5W lPI 



5Ja 



W 



ipi 



but the physical situation is when 

6 a = (h la - (h 2a = 



(33) 



(34) 



Fourth order full quantum correlations from a Langevin-Schwinger-Dyson equation 7 
When this obtains, then 

a + = \ la + 2a ) = r (35) 

is the physical expectation value 

The 1PI effective action is the Legendre transform 

Tip/ [<p A ] = W 1PI - J A <f) A (36) 

Therefore 

S^ipi 



1PLA 



-Ja (37) 



The 1PI CTP EA can be written generically as p] 



T [0_, 0+] = a _D a [0+] + ^f_N ai [0+] b „ + ... (38) 

Where the ellipsis means terms of higher order in 0_. Here D (N) is the so-called 
dissipation (noise) kernel. This is of course identical in form to the effective action for a 
stochastic theory [221 an d thereby we may write down the equivalent Langevin equation 

D a [$.] = -£„ (39) 

where the £ a are Gaussian stochastic sources with zero mean and self-correlation 

[iaii) = N a6 (40) 



2.3. Full Hadamard propagator from the stochastic 1PI approach 

As an application of the stochastic 1PI approach we shall show that the two point 
function for the stochastic theory [27] is identical to half the Hadamard propagator of 
the full quantum theory. 

The full quantum propagators may be derived from the 1PI generating functional 

as 

G AB = -i|^f = -Mtff (41) 
oJ A dJ B 

From the properties of the Legendre transform we have 

Tipi,abW^ = -5 C A (42) 

whereby 

T 1pi ,abG bc = i5 c A (43) 
In the A = (a, a) representation, where a = ±, we find 

b c _ [D C)Q + fN cb , a 



^lPI,(a,a) 



D a [0+] + iN ab <j) 



i, 



(44) 
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And therefore the Hessian, evaluated at the physical point 4> c _ = is 

D 6>a 



rW,(a,a),(/3,b) = ( ^ ..-J" j (45) 



We also identify 



GM , m = , jof , (46) 



ah 



adv 



where G^ b t is the retarded propagator, G ab dv = G b r a et is the advanced propagator, and 
is the Hadamard propagator 

G? = ({$•, $y) (47) 

The equations for the propagators become 

/ D 6>a \ / \G b { -iG& K 



, - z f4^ 

D a , b zN afe ; v -zc^ o ; V 8% ) 1 



namely 



D a ,bG r c et — - 8° a 
Db, a G b ad V = — 8° 

B a , b G\ c = - 2N ab G b ; dv (49) 

This last equation shows that indeed G\ is twice the stochastic propagator G ab defined 
in eq. [271 

2.4- Remarks 

In this Section we have shown that the 1PI EA derived from a quantum field theory, cut 
off at terms quadratic in the difference field </>_, is identical in form to the effective action 
derived from a suitable Langevin equation. This allows us, for example, to compute the 
full Hadamard function from the stochastic approach. 

It would be desirable to extend this equivalence to higher correlations, and to this 
effect we may wish to find stochastic equivalents to the higher effective actions to be 
introduced in next Section. However, the stochastic approach as presented above relies 
heavily on the representation [33 for the 1PI EA. The structure of the 1PI EA is not 
replicated in the higher effective actions (except under restrictive conditions on the 
propagators, see [48]) and so this method fails in general. 

In the following we shall present an alternative derivation of the stochastic approach 
which does not rely in any particular feature of the 1PI EA, and therefore it is readily 
generalized to higher effective actions. As a first step, we shall briefly introduce the 2PI 
effective action. 
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3. The 2PI EA and early stochastic formulations 
3.1. The 2PI EA 

The 1PI generating functional introduced in last Section admits a path integral 
representation in terms of fields defined on the closed time-path 



[Ja] = J D $A j[s[3> A ]+J A 3> A ] 



(50) 



where S 



= S f^ 1 ] — S [$ 2 ]* is the classical closed time-path action. In the 2PI theory 
we add nonlocal sources K AB coupled to (1/2)$^$^. Thus the CTP 2PI generating 
functional is 



e iW[J A ,K AB ] = J D ^A e i[s[<f A }+J A <f A +^K AB <f A <f B } 

The first and second derivatives of the 2PI generating potential W read 

W A = cf) A 



W ,(AB) = 1 

2 



lA a B 



G 



AB 



W 



AB 



%G 



AB 



W ,A{BC) = 1 - 4> A + G BC ) 



*h*h*Z*h) ~ (<t> A <t> B + G 



W ,(AB)(CD) 

Observe that in these equations we are writing 



AB\ ( 1C1D 



G 



CD 



5K 



AB 



to distinguish it from 



W 



AB 



5 2 W 



SJaSJb 

The 2PI CTP EA is the full Legendre transform 



r = w - j A <p A - -K AB + G 

Therefore the mean field equations of motion are 
I\ a = -Ja - K AB( p B 



AB 



5T 



1 



r,(AB) = JG 1 ^ ~ ~2 Kab 
One further variation yields the identities 



,AC 



2T 



(AC) 



lC,E 



1 ,A(CD)G = -d A 



(51) 

(52) 
(53) 
(54) 
(55) 

(56) 
(57) 

(58) 

(59) 

(60) 
(61) 

(62) 
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\ AC - 2 r + r, A(CD) G CD ^ = -5[ E / B y 

r,(AS)C0 C ''' B + ^,(AB)(CD)G CD ' E = 

r AC,(EF) , r n CD,(EF) _ 1 

where 5(ab) stands for the symmetrized identity operator 

We now write the derivatives of the mean fields in terms of correlations 

<f>C,E = W ,CE = lQ CE 

^,(EF) = W ,C(EF) = 1 ^C^E^ _ C ^ E(j) F + Q EF^ 



(63) 
(64) 

(65) 

(66) 

(67) 
(68) 



QCD,E = 2W ,(CD)E _ ,C,E,D _ m^C 



^>l) ~ <P E <P C <P D ' <P E G CD - <f> C G DE - <P D G CE 



(69) 



qCD,{EF) = 2W ,(CD)(EF) _ ^CiEF)^ _ ,(E F) 

= \ - (<P C <P° + G CD ) {f<f + G EF ) 

-<j> D [(^K)-^(^V + g bf )" 

This last equation may be rewritten as 



(70) 



_ ^ (j) C (j) D + G CD^ (j) E (j) F + G EF^ 



and further 



(71) 



- 2tG CD ^ = 

- <f<fP<f> E <f 

- <f<fG EF - <p c <p E G DF - <p c <p F G DE 

- <t> D <t> E G CF - <p D <p F G CE - <t> E 4> F G CD 
+ i<j> c G EF ' D + i<l> D G EF ' c 

We also notice the identity 

- 2i<p c ^ = -iG CE > F + <P E G CF + <P F G CE 



(72) 



(73) 
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3. 2. Early stochastic formulations 

As we have said in the Introduction, it has been known for a long time that the 
Boltzmann equation is just a mean field equation, and can be improved by upgrading 
it to a full Langevin type equation where particle number fluctuations are explicitly 
included. Since from the point of view of field theory the Boltzmann equation is 
just a particular limit of the Kadanoff-Baym equations, which are in turn equivalent 
to the Schwinger-Dyson equations [1], it is natural to seek a corresponding stochastic 
generalization of the latter. In this subsection we shall review two early attempts in this 
direction and point out their shortcomings. The proper stochastic formulation shall be 
presented in next Section. 

As we have said in Section II, the peculiar structure of the 1PI EA allows to associate 
to a field theory problem an equivalent Langevin equation, whereby, for example, the 
Hadamard propagator may be obtained as a stochastic average. One possible strategy is 
to try to cast the 2PI EA in a similar framework. This line of thought is pursued in [48J. 
Success is found only under special, and restrictive, assumptions on the structure of the 
propagators, and therefore it is unsuitable as a general foundation for the formalism. 

In [49] the same authors follow a different strategy, which may be regarded as a 
nonequilibrium generalization of the fluctuation - dissipation theorem. We shall consider 
here only the case of a field theory with no background mean fields. The general case 
will be discussed in next Section. 

We assume there are stochastic kernels G AB = G AB + 7 such that the stochastic 
averages of these kernels reproduce the quantum averages of products of the composite 
operator These kernels obey the Langevin equation 



We have two ways of computing the self correlation for the stochastic kernels. By 
assumption 



T (AB) [G s ] = —KAB 

After linearization, this becomes 



(74) 



r (AB)(OD)7 = ~7T K AB 



(75) 



(G AB G C S D ) = (^>h^h) = G AB G CD - UW'WW 



(76) 



while from the explicit solution of the linearized equations we get 



(G AB G C S D ) = G AB G CD + { 1 AB 1 CD ) 



= G G + - T t (AB)(EF) F,(CD)(GH) (^EF^Gh) (77) 

Asking both computations to agree we get 



^i^,(IJ)(AB)^,(KL)(CD) 



W ,(AB)(CD) = 



,(IJ)(KL) 



(78) 
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As we shall see in next Section, this is the correct result. However, the authors of 
provide some contrived argument to the effect that the sign of the last term in eq. (1751) 
ought to be changed. These arguments would be unconvincing, except the sign reversal 
seems to be necessary to derive the correct noise term for the Boltzmann - Langevin 
equation in the kinetic theory limit. 

This seeming paradox demands clarification, which will be provided in next Section. 

4. Systematic Langevin approach 

4-1- Systematic Langevin approach to mean field dynamics 

So far in this note we have seen two different strategies to associate an equivalent 
stochastic equation to a field theory problem. The approach of Section 2 depends upon 
the formal analogy between the 1PI EA for the quantum and stochastic problems. 
This approach is straightforward for the 1PI theory, but it is not easily generalizable 
to higher effective actions. The approach of Section 3 assumes there is an equivalent 
stochastic description, where the noise is whatever it needs to be to sustain the proper 
fluctuations. In this sense, it is a generalized fluctuation - dissipation kind of argument. 
This approach works for any effective action, but seems to be subject to a sign ambiguity. 
In this Section we shall develop this second approach in a systematic way, dispelling the 
apparent paradox afflicting the discussion in |49j. 

Let us reformulate the 1PI problem in terms of the fluctuation - dissipation inspired 
approach. To derive the Langevin-like approach, we write down eq. (J371) for a stochastic 
c-number CTP field variable $^ coupled to a stochastic source £^ p/ . We assume the 
solution takes the form 



where the displacements <p A account for both the quantum and statistical fluctuations 



(79) 





In the linearized regime, we get 



^lPI,AB<f 



B 




We find 





(83) 



so we require, using that G = G 



'DB 





that is 
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It is important to stress the differences between this approach and that of Section 2. 
In Section 2 there was only one stochastic source coupled to the stochastic condensate 
field. In the new approach of this section, the number of stochastic sources has been 
doubled. 

Other way to express the same concept is to observe that the constraint 0_ = 
which characterizes physical mean field configurations is violated in our stochastic 
formulation. The formal stochastic equations for the fluctuations read 



V D b,a = -£+a 

D a ,^ +6 + *N afcV ^= -e P J (86) 
The first equation yields 

^ h = G b : dv e + 7 (s?) 

We can now eliminate (p~ b and obtain a single stochastic equation for ip +b . The question 
arises if this single equation is the same as derived in Section 2. This requires we identify 

I = £7 + MabGlU+c 1 (88) 

Now we have = -N ab , (£™&£ J ) = *D a , b and ) = 0, so the above 

identification implies 

(Lib) = -N„5 - N M GlDa, c - N fld GlD 6 , c = N ab (89) 

as expected. After the auxiliary field ip~ b has been integrated out, both stochastic 
formulations are equivalent. 



4-2. The stochastic 2PI theory 

One advantage of the new formalism is that it can be trivially generalized to higher 
effective actions. Let us consider again the 2PI theory. 

This time we shall consider the general case where symmetry may be broken. To this 
end we write as before <3>^ = (f) A + ip A , but now we also write a stochastic representation 
for the composite operator — (f) A ^j — 4> B ^j . In the stochastic theory this 

composite operator is represented by a stochastic c- number kernel G AB = G AB + ^ AB . 
We assume the identity 

$^$f = + G AB + + V A <P B + 1 AB (90) 
The required expectation values are 

= <f> C <f> E + <>V) (91) 



= tpC^D^E + ^CqDE + ^DqCE + ^EqCD + ^DE^ (Q2) 
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Observe that this implies that (^ c ^i DE ) is totally symmetric. 



= + G CD ) (0 V + G 



+ 



FQDF + qDE + ^^ 7 ^y 

+ £ (^ F 7 CD ) + F (^7 CD ) + {l CD l EF ) 



(93) 

We can now relate the derivatives of the mean fields with respect to the sources to 
stochastic averages 



!) C > E = i (ip C ip E 



qCD,e 



C„,DE\ 



<P C ^ = - {( V c j EF ) + <P E G CF + <p F G CE ) 



(94) 
(95) 

(96) 



qCD,(EF) = « | 0£ ^Jy^ + f 



^V D ) + (7 CD 7 £F >} 



Whereby we get the identities 



/V) + r A(CD) ( 7 ^> 



,CD,„E 



i5 E 



(97) 



(98) 



I> - 2T {AO) ] {(^ C 1 EF ) + <P E G CF + <P F G CE 



+ r A(CD) uz/,^,co\ 



which reduces to 

T AC — ^,(AC) 



v f i cd ) + F (<p e i cd ) + (i cd i ef )} = 2<jy 



,-C-,EF\ , r /„,CDEF\ n 

V 5 7 / + 1 ,A(CD) \7 7 / = 



r )( AB)c V) + r )( AB)(CD) (7 Ci V) = o 



(99) 

(100) 
(101) 

(102) 



r /,„C n ,EF\ , t / ',CD n ,EF\ 1 / rBrF , xFcE\ 

r ,(4B)C^ 7 / + r,(AB)(C£») \7 7 /=2l AB+ AB J 

Assuming the mean field equation F,(ac) = these equation suggest a stochastic 
dynamics for the ip, 7 fields 

(103) 



r Acy? c + r a(cd)7 cd = -£4 



Observe that a possible term kab(P B is absent 



T (A^CV^ + r j (AB)(CD)7 C ' D = ~q^ Kab 



(104) 
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provided 

[uv E ) = -«5 ( 105 ) 

[Uf°) = (106) 

'K AB ip E ) = (107) 

[kab1 EF ) = ~i {Sffi + (108) 

Multiplying the Langevin equations by the sources and using these expectation values, 
we get 

(USb) = ^T,AB (109) 

(Za^cd) = 2iT A{CD ) (110) 

(KCDKAB) = MY ,{AB){CD) (HI) 

4-3. Recovery of the 1PI stochastic theory from the 2PI one 

Let us check that the 2PI and 1PI theories agree as far as the mean field fluctuations 
are concerned. This must be true, as the 2PI theory is built around the requirement 
that (ip A y B ) = G AB 

The 1PIEA Tipi is recovered from the 2PIEA T as 

r 1PJ [0] = r[0,G o [0]] (ii2) 

where the correlations Go are slaved to the mean field through 

r (AB )[0,G o [0]] =0 (113) 
One further derivative shows that 



I", (AB)(CD) 



-1 



f,(ab)e (114) 

Therefore the first and second derivatives of the 1PIEA are 

T lPI , A [0] = T )A (115) 

r iP i,AB [<P] = T, AB + T,A(CD)G^ (116) 

The stochastic equation for the mean field fluctuations, as derived from the 1PIEA, is 
simply eq. flHTj) with noise self-correlation fl85l) . Now eq. fjl04j) admits the solution 

„,CD n CD,E 1 Pp ]' 1 ^ (M7\ 

7 = G Q E Lp - - [1 ,(AB)(CD)\ I^AB (117) 

Using this into eq. (11031) we get back eq. (J8T]) provided we identify 

e A PI = U + \g c ^cd (118) 
Indeed, both expressions have the same self-correlation. 
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4-4- Stochastic equations for the physical propagators 

Finally, let us recover the noise self-correlation given in [49]. The situation is similar 
to the one in the 1PI theory. In the 1PI theory, the <p~ field vanishes identically on- 
shell. However, in the stochastic approach we assign a nontrivial source £i P/ to it. It 
is by eliminating this auxiliary field that we recover the usual approach, with a single 
stochastic source £ whose self-correlation is given by the noise kernel. 

Similarly, in the quantum field theory problem the correlator G = (ip~(p~) 
vanishes identically, as a result of path ordering. However, in the stochastic approach, 
we consider it as an auxiliary field and couple a source to it. The authors of |49j failed 
to recognize the violation of the constraint G = 0, but compensated this oversight by 
forcing a sign change in their expression for the noise self-correlation. In the context of 
the proper theory, this sign change is due to the elimination of the auxiliary field G , 
keeping only the physical degrees of freedom. 

For simplicity, and as in [49J , let us assume a symmetric scalar field theory, meaning 
that the background fields vanish and also T^bc) = 0- Then the 2PI EA can be written 
as 



2 — 2 , (119) 

The first derivatives of the 2PI CTP EA yield the mean field equations of motion 

S AB - iG~ AB + 2V Q>(AB) = (120) 

we shall call G AB the on-shell propagators. We adopt the (±, a) indexes, as before. 
In this representation, we have the identifications eq. fj4*6|) and (jUD - We now replace 
the generic kernels G AB by the stochastic kernels G S AB = G AB + 7 AB and expand the 
effective action to second order. For simplicity, we only keep linear terms in Tq. We 
show in [49J that this is enough to study the fluctuation terms in the kinetic field theory 
limit. The relevant quadratic terms in the effective action are 



r (2) 

where, as in eq. 

G 



4 



-Tr 





therefore 



and 



-*D M7 (- fe )(+ c ) 

-*D a , b7 (+6){+c) + N ab7 M){+c) 



-*D a , fe7 (+b)( ~ c) 



(121) 

(122) 

<-&)(-c) \ 



(123) 



r (2) = ^ {d m7 (-»(+ c )d 



4 

+ 2D b>o7 



(-d)(+a) 



(-b)(-c) 



(+d)(+a) 



+ iN cd7 



(-d)(+a)" 



+ 



D a,fe7 



(+b)(-c) 



(-6)(-c)l 



D c,d7 



(+d)(-a) 



(-cO(-a)" 



)(124) 
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The stochastic equations for the propagators are 



r 



^ +c >D dc + n ( - b)( - c) N 



cd 



■«(-«0(+a) 



—I 



D c ,d = — K 



(+d)(-a) 



(125) 
(126) 



^ {D M [D c , d7 ( +d )( +a ) + iN*^"^ 



2 

+ ^N a6 



-1 



B c , d ^(-») +zN cd7 (- d )(-°)]} = -^K(_6,(_c) 



This last equation implies 

(-fc)(-c) _ _;ncd nba 
1 — l(j adv (j adv K {+d){+a) 



(127) 
(128) 

(129) 



Observe that not zero in the equivalent stochastic problem. The equations 

for the physical stochastic propagators are obtained eliminating 7 (~ 6 X~ C ) throughout. 
For example, for the fluctuations in the Hadamard propagator we obtain 



D c , d7 ( +d )(+ a )D M +zN cd7 (- d )( +a )D M +zD Ci(i7 ( +d )(- a )N afe = -m ( _ 6)( _ c) (130) 



where 



«(-6)(-c) = K(-6)(-c) +^N cd7 (- d )(- a )N ab 

= K(-6)(-c) + N ceG^G^ t ,N /fe K( +d) ( +a ) 

From the self-correlations of the original sources 

K(-6)(-c)«(-d)(-a)) = — 2N a6 N C(i 
K (-6)(-c)^(+d)(+a)) = 2D bi(l D C:d 
K (+6)(+c)^(+d)(+a)) = 



we get 



(_6)(-c)«(-c0(-a) 



2N..N 



ab^cd 



(131) 

(132) 
(133) 
(134) 

(135) 



which agrees with |39] . Observe that in going from eq. f 11 3 2 j) to eq. (11351) the sign of the 
right hand side has changed. Thus we obtain the sign change enforced in [49J, whose 
origin is now clear. 
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5. Final remarks 

In this paper we have presented a systematic stochastic approach which allows to 
build Langevin-like generalizations of the Schwinger-Dyson equations of an interacting 
quantum field theory. The resulting fluctuations account both for quantum and 
statistical fluctuations. 

It is hoped that the present approach will prove clearer and easier to apply 
than existing alternatives in the literature. Quantum and statistical fluctuations 
are increasingly under scrutiny in fields that range from cosmology to Bose-Einstein 
condensates. We expect the tool we have presented in this note will prove valuable in 
this endeavor. 
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